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FEASIBILITY OF SUPPORTING LIQUID FUEL ON A SOLID WALL 
IN A RADIATING LIQUID-CORE NUCLEAR ROCKET CONCEPT 
by Henry A. Putre and Albert F. Kascak 
Lewis Research Center 

SUMMARY 

A reference liquid fuel nuclear rocket concept, previously analyzed, consists of a 
liquid carbide fuel mixture supported on the inside wall of a rotating tube. Most of the 
nuclear heat generated in the liquid fuel is radiated to the axially flowing hydrogen pro- 
pellant. A fraction of this heat generation, called the wall cooling ratio, is removed by 
support-wall cooling. The reference engine delivered a 1430-second specific impulse, 
with a 9500° R (5280° K) radiating surface temperature, and a 25-percent wall cooling 
ratio. Solid-support-wall temperature for this rocket, assuming only simple molecular 
conduction in the liquid fuel, was estimated as 14 500° R (8000° K), which of course, 
greatly exceeds structural material limits. 

This report extends the original liquid fuel heat-transfer analysis for the reference 
engine by developing an expression for combined molecular and propellant induced turbu- 
lent heat transfer in the liquid fuel, and results in a more rigorous determination of 
support-wall temperature. In addition, two means of reducing support -wall temperature 
without severely affecting rocket performance are investigated. 

For the reference fuel-element configuration, with uniform heat source distribution 
in the liquid fuel, the support-wall temperature calculated from the more complete tur- 
bulent analysis was 9830° R (5460° K), which is well above currently practical wall ma- 
terial limits. 

When the fuel was zoned into a region in the liquid near the propellant interface, the 
support-wall temperature was found to be about 1200° R (670° K) lower than for the un- 
zoned case. A simplified mass-transfer analysis indicates, however, that fuel zoning 
cannot be maintained over rocket lifetime, and therefore wall temperature reduction by 
fuel zoning does not appear feasible. 

The insulating effect of a stable opaque vapor layer built up at the solid-liquid inter- 
face was also evaluated. Wall temperature was found to be reduced to 8500° R (4720° K) 
at 1430 seconds specific impulse by the vapor layer. Current wall materials, such as 
graphite used at 7200° R (4000° K) in conjunction with an opaque vapor layer, may per- 
mit operation at a reduced specific impulse of about 1250 seconds. 



INTRODUCTION 


For planetary space missions and short trip times, a high specific impulse and mod- 
erate thrust rocket is necessary. Reference 1 describes such a nuclear rocket, which 
employs a bundle of rotating liquid fuel elements with the liquid niobium carbide (NbC)and 
uranium carbide (UCg) mixture supported on solid cylindrical walls by centripetal accel- 
eration. Figure 1(a) shows one of the fuel elements. The propellant flows axially through 
the center of the annular liquid fuel element, and it is heated by thermal radiation from 
the liquid surface. 

The ultimate goal of this rocket is to heat the propellant to as high a temperature as 
possible with minimum mass transfer from the liquid fuel element to the propellant. The 
fuel temperature is chosen sufficiently high so that thermal radiation is the dominant 
mode of heat transfer, thus almost uncoupling heat and mass transfer. 

Reference 1 showed that from cycle considerations, the percent of heat removed 
through the solid wall by regenerative cooling strongly affected the rocket specific im- 
pulse. In fact, for maximum specific impulse this wall cooling ratio must be kept to a 
minimum (see table I). With reduced wall cooling, however, the temperature of the 

TABLE I. - WALL COOLING RATIO 
AND SPECIFIC IMPULSE FOR 
3500° R PROPELLANT 
TEMPERATURE AT 
FUEL-ELEMENT 
INLET (FROM 
REF. 1) 


Wall cooling 
ratio, 

V 

Specific impulse, 

I , 

sp’ 

sec 

0.2 

1570 

.25 

1430 

.3 

1280 

.4 

1110 

. 5 

1000 


inner solid wall surface tends to increase above the wall material structural limit. The 
feasibility of this rocket concept then rests on the possibility of finding a compromise 
between these two limits. 

The reference engine evolved in reference 1 had a 1430-second specific impulse with 
a radiating surface temperature of 9500° R (5280° K). Twenty -five percent of the nuclear 
heat generation had to be removed by cooling the solid wall and surrounding structure. 
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An estimate of solid support-wall temperatures for this engine, considering only molec- 
ular conduction in the liquid fuel, indicated that the wall temperature (14 400° R, 8000° K) 
greatly exceeded structural material limits and that the fuel would boil internally. 

The analysis in reference 1 had several shortcomings, but it did serve to show that 
turbulent effects in the liquid fuel mixture should be considered in a determination of 
support -wall temperatures. The objectives of this report are to extend the liquid fuel 
heat-transfer analysis of reference 1 by including turbulence, and to determine accu- 
rately support-wall temperature at design performance conditions. Specific questions to 
be answered in this report are the following: 

(1) With turbulent convection in the liquid fuel properly accounted for, do support- 
wall temperatures fall within realistic material limits? 

(2) Are support -wall temperatures significantly reduced by zoning the fuel into a re- 
gion in the liquid away from the support wall? 

(3) Could an opaque vapor layer built up at the liquid-solid interface reduce support- 
wall temperature to within realistic limits? The maximum stable vapor -layer thickness 
is estimated from a simple vapor bubble force balance. 

With these objectives in mind, the liquid-fuel -element configuration of reference 1, 
illustrated in figure 1(b), was analyzed. Design values from reference 1 which will be 
used in this study are given in table II. 


TABLE n. - DESIGN VALUES FOR LIQUID-CORE 
NUCLEAR ROCKET FROM REFERENCE 1 


Fuel mixture temperature at gas surface, 

9500° R (5280° K) 

Radiant heat flux from liquid -fuel mixture to gas stream, q” 

3. 52xl0 7 J/(m 2 )(sec) 

Rocket specific impulse, I 

1430 sec 

Wall cooling ratio, rj 

0.25 

Thickness of liquid -fuel mixture, L 

0. 0095 m 

Inner radius of fuel -element solid support wall, R 

0. 060 m 

Rotational speed of fuel -element support wall, w 

2400 rpm 

Radial acceleration due to element rotation, g r 

378 g Q 

Propellant mass velocity, 

33. 0 kg/(m 2 )(sec) 

Propellant density, 

1. 47 kg/m 3 


SYMBOLS 

A ratio of eddy to molecular thermal conductivity at y = L 
a^ Rosseland mean radiative absorption coefficient 
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c 

CP L 

D 

d 

f g 

°g 

So 

gf 

*sp 

K(r) 

k L 

L 

Pr L 

Q 

q”’(r) 

*1 

R 

r 

S 

T 

T 


w 


1/2 


mole fraction of uranium carbide in liquid fuel mixture 
value of C at solid wall 

value of C averaged over the liquid fuel thickness 
liquid specific heat 

molecular mass diffusion coefficient 

diameter of vapor bubble at solid wall, see fig. 3(a) 

vapor bubble base diameter, see fig. 3(a) 

pipe friction factor for gas stream 

gas mass velocity 

2 

gravitational acceleration constant, 9. 80 m/ sec 
radial acceleration 

specific impulse, defined as thrust per unit weight flow rate at g Q , sec 

generalized total thermal conductivity 

molecular thermal conductivity 

thickness of liquid fuel mixture 

liquid Prandtl number 

constant for nuclear heat generation rate per unit volume 
heat flux from liquid fuel mixture to solid wall 
generalized nuclear heat generation rate per unit volume 
heat flux from liquid fuel mixture to gas stream 
inner radius of fuel -element solid wall 
distance from fuel -element centerline 
equivalent surface roughness 
liquid -fuel -mixture temperature 

fuel -mixture surface temperature facing the solid wall 
fuel-mixture temperature at gas surface 
solid-wall temperature at inner surface 
time 

real time corresponding to $ 1/2 
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liquid-fuel -mixture friction velocity 


distance in liquid fuel mixture measured from solid wall 

y + = (yu tfj VjJ, dimensionless value corresponding to y 

dimensionless pipe roughness parameter, a = S/L 

fraction of liquid thickness which has no heat sources 

vapor -layer thickness 

liquid eddy thermal diffusivity 

liquid eddy momentum diffusivity 

ratio of heat removed at solid wall to total heat generated 

9 

dimensionless time, 6 = Dt/L 

dimensionless time when C /C = 0. 5 

o 

liquid viscosity 

liquid kinematic viscosity, p^ = Mj_/Pl 
gas density 
liquid density 

liquid -fuel -mixture surface tension 
gas shear stress at liquid -gas interface 
contact angle at base of vapor bubble, see fig. 3(a) 
optical thickness, number of radiation mean free paths 
fuel-element angular velocity 


ANALYSIS 

Since this report treats several separate but related problems, as indicated in the 
INTRODUCTION, the analysis of these problems is presented in several separate sec- 
tions. The relevancy of each of the various analysis sections is discussed in greater 
detail as part of the RESULTS AND DISCUSSION. 

The first analysis derives the temperature distribution in the liquid fuel mixture, 
with turbulent heat transport included, for a generalized step distribution of heat sources 
due to nuclear fuel zoning. The next section (p. 17) derives the fuel concentration distri- 
bution in the liquid mixture as a function of time from rocket startup for various cases 
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of initial fuel zoning. The final section (p. 18) uses a radiative heat-transfer analysis to 
calculate the temperature drop across a stable, insulating, opaque vapor layer built up 
at the liquid-solid interface. 


LIQUID-FUEL-MIXTURE TEMPERATURE PROFILES WITH ZONED 
HEAT SOURCES AND TURBULENT HEAT TRANSPORT 

The temperature profiles in the liquid fuel mixture and the wall temperature T w are 
calculated for the fuel element model illustrated in figure 1(b) using the reference design 
radiating surface temperature and heat flux and the range of wall cooling ratios in refer- 
ence 1. This calculation requires that heat source and thermal conductivity distributions 
are specified throughout the liquid. The heat diffusion equation with only radial steady - 
state heat transfer is considered. In the most general case, both the heat source strength 
and thermal conductivity yary radially. The generalized thermal conductivity includes 
both molecular and eddy heat transport. 

The generalized heat-transfer relation for the liquid fuel mixture then is 

— — frK(r) — 1 + q M? (r) = 0 (1) 

r dr L drj 

The boundary conditions from figure 1(b) at r = (R - L) are 

T = Tj = 9500° R (5280° K) (2) 

K(r) — = q’’ = 3. 52xl0 7 (3) 

dr ^ p 

ar (nO(sec) 

(These design values are given in table II. ) 

Temperature profile solutions are seen to be physically realistic to the extent that 
expressions for q M, (r) and K(r) apply. Thus these quantities must be carefully con- 
sidered. The heat source term q"'(r) is by definition related to the wall cooling ratio 
tj and the gas interface heat flux qy by the equation 



R 

- L 


rq"'(r) dr = q" (R - L) + q"R = qV 
1 oi 



(4) 


For the case of zoned nuclear fuel, a step heat source distribution is assumed with a 
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uniform nonzero heat source in the region /3L < y < L, as indicated in figure 1(b). With 
this zoned heat source distribution and equation (4), expression for equation (1), rewritten 
in terms of 77 and qy with y measuredfrom the support wall, becomes, for 0 < y < /3L, 

q"»(r) = 0 (5) 


and for |3L < y < L, 


q'"(r) = Q = 

(! - v) 


2q’ 1 ’(R - L) 

_ c •- - ^ • 

(R - /3L) 2 - (R - L) 2 


( 6 ) 


An important phenomenon now considered is that of turbulent mixing in the liquid 
fuel mixture. The objective is to derive expressions for shear stress and finally eddy 
diffusivity in the liquid fuel mixture. From studies of liquid-gas surface interaction in 
annular two -phase flow, it is well known that, with a turbulent gas stream flowing over a 
thin liquid film, the turbulent gas eddies interact with the liquid surface. They thereby 
transfer their random mixing motion to the liquid. By postulating that the liquid-film 
surface in annular two -phase flow affects the gas like a rough -wall pipe, Wrobel and 
McManus (ref. 2) were able to show that a shear stress could be calculated at the liquid- 
gas interface from known values of gas Reynolds number and liquid -film thickness. Ref- 
erence 2 also showed that this shear stress could be used to find the velocity distribution 
and eddy viscosity in the liquid film. 

Similar liquid-gas surface interactions are assumed to occur in the liquid -fuel - 
element configuration; then, turbulent mixing in the liquid fuel mixture is an important 
mechanism for heat transfer. An expression for liquid fuel eddy thermal conductivity is 
derived with the aid of the analysis of reference 2. 

The liquid fuel acts on the propellant gas like a rough-wall pipe having an equivalent 
roughness equal to some fraction a of the liquid thickness, or 

S = ah (7) 

With this surface roughness used in a rough -wall -pipe friction factor correlation and 
with the hydrogen propellant velocity given in table II, the interfacial shear stress can be 
calculated. The friction factor correlation for rough-wall pipes is given in reference 3 
(p. 525); for the gas stream it takes the form 


7 



1 


( 8 ) 


2 log ( R - L j + 1.74 

The correct value of a to be used in equation (7) depends on the wave structure of 
the liquid surface. Based on the work of reference 2, a value of a = 1/2 gives at worst 
an underestimate of gas friction factor and liquid turbulence level, and therefore, it is 
used here. A more accurate value of a requires detailed experimental study of the gas- 
liquid surface interaction. Such a study was not considered here since the end results 
are quite insensitive to the exact value of a. 

The gas friction factor calculated from equations (7) and (8), using design values 
from table II, becomes 


f g - 0. 070 (9) 

The gas shear at the gas -liquid interface in terms of gas mass velocity and friction 
factor is 


T 


f G 2 

J g. 

8p 

S 


( 10 ) 


The liquid fuel eddy viscosity can now be calculated assuming this interfacial shear 
value applies for the entire liquid thickness, and that this shear stress leads to turbulent 
mixing in the liquid fuel similar to mixing in developed pipe flow for the same shear 
stress. The eddy viscosity can be derived from the universal velocity distribution given 
in reference 3 (p. 509) with the following results: 

For y + < 5 (laminar sublayer), 


e 


m 


= 0 


( 11 ) 


For y + S: 70 (turbulent region), 


e 


m 


= ^ L (0.4y + - 1) 


( 12 ) 


where 
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(13) 


y 


+ 



and 



(14) 


For simplicity and for purposes here, equation (12) is assumed to apply over the re- 
gion 5 < y + < 70. Invoking the Reynolds analogy between turbulent momentum and heat 
transfer, or = e , gives the expression for thermal conductivity including turbulent 
mixing as 


K(r) = k L y + < 5 (15) 

K(r) = k L [l + Pr L (0. 4 y + - 1)| y + > 5 (16) 

Using liquid-fuel-mixture properties from the appendix and design values from 
table II enables equations (15) and (16) to be rewritten in the nearly equivalent form 

+A (S) (17) 

where the calculated value of the turbulence coefficient is 

A = 11.0 (18) 

This variation of total thermal conductivity across the liquid thickness is shown in fig- 
ure 2. As is generally found in turbulent boundary layer flows, total thermal conductivity 
increases linearly from the solid wall (see fig. 2). 

Using the total thermal conductivity in the form of equation (17) together with the 
heat source terms of equations (5) and (6) in the basic heat-transfer equations (1), (2), 
and (3), and satisfying temperature and heat flux continuity at y = j3L yield the tempera- 
ture solution for the zoned fuel case with turbulent conduction: 


K(r) = k L 1 
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For 0 < y < /3L: 


T = T(y = /3L) - 


QrjL 

2k T 


(R - /3L) 2 - (R - L) 2 


L + AR 


In 


1 + Al' 
R-y L 


- In 


R - /3L 1 + Aj3 


(19) 


For /3L < y < L: 



9 9 

1 + A L 

QL 

(1 - V )(R - /3L) 2 + tj(R - L r 

In R - y - m L 

2k L 

L + AR 

R - L 1 + A 



— 


A'1“l) + (i + l“ 


uaI 

AR\, L 


1 + A 


( 20 ) 


where 


T(y = /3L) = Tj_ + 


QL 

2k T 


(1 - 77 ) (R - £L) 2 + 77(R - L ) 2 


L + AR 


ln R -^_ ln I±Ag 

R - L 1 + A 


A(1 - 0) + (l + ^ln 1 + . 

\ L / 1 + A 


( 21 ) 


TIME DEPENDENCE OF ZONED FUEL CONCENTRATIONS 

The time dependence of zoned fuel concentrations is calculated by assuming a simple 
one-dimensional diffusion model. Referring to figure 1(b), and assuming a constant 
value of diffusion coefficient and a small molar fuel concentration compared to sol- 
vent concentration, yields the mass diffusion equation with initial and boundary con- 
ditions for the fuel species as 


d 2 C _ 1 9C (22) 

ay 2 D 9t 


Initial conditions: 
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At t = 0, 


0 < y < 0L (23) 

j3L < y < L (24) 

Boundary conditions: 

At y = 0, 

— =0 (25) 

9y 

and at y = L, 

— = 0 (26) 

9y 

The quantity C in equation (24) is the equilibrium fuel concentration averaged over 
the liquid mixture thickness. The actual value of C is not important here; however, for 
purposes herein, it remains constant for all values of /3. When the separation of vari- 
ables method is used, the solution to equations (22) to (26) is 



£= 1 - 



2 2 

sin(n7r/3)cos^n77 y ' ■'~ 71 71 @ 


. 2 5 

Y_\ e -V * 


n 


where 6 is the dimensionless time variable 


(27) 



(28) 


Equation (27) evaluated at the wall becomes 
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oo 


c 


o 


c 



sin(n7rj3)e -n 71 ® 
n 


(29) 


VAPOR INSULATOR BETWEEN LIQUID AND SOLID REGIONS 

One method of reducing support-wall temperature is by means of an insulating vapor 
layer built up between the support wall and the liquid fuel. Conceptually the vapor layer 
would form either by sublimation of the solid wall or by gas injection through the solid 
wall. The stable thickness of such a vapor layer is difficult to estimate in view of the 
complicated two-phase flow structure; however, an approximate value of the maximum 
stable vapor layer thickness is derived from a single vapor bubble force balance. 

The principal forces acting on a vapor bubble attached to the solid wall in the cen- 
trepetal gravitational field are assumed to be the bouyancy and surface tension forces. 
Then, with a spherical bubble shape assumed as shown in figure 3, the force balance at 
detachment is 


P L — g r = Jra sin (p (30) 

The quantity d^ sin cp is approximated by the bubble diameter. The maximum stable 
vapor thickness is assumed to be equal to the bubble diameter at detachment. 

The maximum stable vapor thickness is thus given by 


6 = d = 



(31) 


Using the value of g r = 387 g Q from reference 1 and the fuel mixture properties given in 
the appendix yield 


6 = 0. 57 mm (32) 

This estimate of the maximum vapor thickness is probably an overestimate since fluid 
dynamical forces on the bubble, which tend to sweep the bubble off the solid wall, have 
been neglected. 


12 


The temperature drop across this vapor layer can be calculated. Since radiation is 
the dominant mode of heat transfer from the liquid fuel to the propellant, only radiative 
heat transfer is assumed through the vapor layer. (This is equivalent to assuming that 
the thermal conductivity of the vapor times the temperature gradient is small compared 
to the amount of heat transferred by radiation. ) The diffusion analysis of reference 4 
for radiative heat transfer across an absorbing medium between flat plates is used, with 
the walls assumed black. This analysis requires that the optical thickness be specified. 

The optical thickness x of the stable vapor layer just described is now estimated 
from the definition 

X 3 J* a R dy (33) 

The determination of x requires that the vapor composition and the absorption coefficient 
of its various species be known. The composition of the vapor is very difficult to deter- 
mine since the diffusion process is nonisothermal. If this vapor has a constant absorp- 
tion coefficient, then x is given by the product of the absorption coefficient and the vapor 
thickness: 


X=a R 5 

1 7 2 

Using the value of a R (3. 2 mm" ) for carbon vapor at 200 atmospheres (2. 02x10 N/m ) 

and 9000° R (5000° K) based on reference 5 gives x = 1* 8. 

It is difficult to assess the accuracy of this estimate of the vapor opacity. Photoion- 
ization and atomic line absorption are normally considered unimportant at temperatures 
less then 10 000° K; molecular transition absorption should be dominant. The absorption 
coefficient is thus characterized by a great number of lines. This creates the possibility 
of a window in the absorption spectrum. The possibility of this window plus the strong 
temperature dependent properties would seem to overshadow the neglect of some of the 
molecular species present; thus, the present estimated x is probably too high. 


RESULTS AND DISCUSSION 

The preceding analysis was applied to the liquid-core nuclear rocket described in 
reference 1, which consists of a liquid carbide fuel mixture supported on the inside wall 
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of a rotating tube. The axially flowing hydrogen propellant is heated by thermal radiation 
from the liquid surface. For purposed here, engine design values were fixed at those of 
reference 1 and shown in table II with the exception of wall cooling ratio and specific im- 
pulse. The wall cooling ratio was varied from 77 = 0 to 0. 5, with primary emphasis 
being on the design value 77 = 0. 25. 

In the following sections the total thermal conductivity, given by equation (17), was 
used to evaluate liquid -fuel -mixture temperatures. The first section discusses the im- 
portance of turbulence and evaluates these temperatures for the case of uniform heat 
source distribution. The second section discusses the feasibility of fuel zoning, and the 
last section discusses the advantages of an insulating vapor layer between the liquid fuel 
and the solid support wall. 


Fuel-Mixture Temperatures With Uniform Heat Source Distribution 

Temperature profiles in the liquid fuel mixture for uniform heat source distribution 
and assuming only molecular conduction are shown in figure 4. These profiles were cal- 
culated using equation (20) with /3 = 0 and A = 0. Two important inconsistencies occur 
when only molecular conduction is considered: 

(1) Since the boiling point of niobium carbide, the principal fuel -mixture component, 
is approximately 13 000° R (7200° K) (see the appendix), figure 4 indicates that fuel boil- 
ing occurs in the liquid, resulting in an additional heat-transfer mode. 

(2) The large radial temperature gradients shown in figure 4 coupled with a large 
radial acceleration field due to fuel-element rotation will cause free convection effects in 
the liquid near the propellant interface. This results in yet another additional heat- 
transfer mode. Both of these additional heat-transfer modes are ignored by an analysis 
based only on simple molecular conduction. 

Temperature profiles for the same case but with turbulence included are shown in 
figure 5 and were calculated using equation (20) with /3 = 0 and with the calculated value 
A = 11. 0. A comparison of figures 4 and 5 indicates that liquid fuel temperatures near 
the support wall may be greatly overestimated when turbulence is ignored. The eddy 
conductivity is generally larger than the molecular conductivity (see fig. 2), so that both 
the temperature gradient and the maximum temperature in the liquid are significantly 
reduced. Thus, in the presence of eddy conduction, internal boiling will probably not 
occur. Also, the importance of free convection is greatly reduced. Therefore, when 
turbulence is accounted for, the two inconsistencies present in the laminar calculation 
are much less pronounced, or disappear altogether. 

Support -wall temperatures for the uniform heat source distribution and various wall 
cooling ratios were calculated using equation (20) with /3 = 0, and they are illustrated in 
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figure 6. Values of wall temperature are plotted for the calculated turbulence coefficient 
A = 11.0 and also for the A = 5. 0 and A = 20. 0. Figure 6 indicates that, using the tur- 
bulence theory, wall temperatures are quite insensitive to the actual value of A. Thus, 
any uncertainties in fluid properties, or applicability of the turbulence theory, which may 
cause a deviation in the turbulence coefficient of as much as a factor of 2 from the calcu- 
lated value A = 11.0, will not greatly affect the nature of the wall temperature curve. 

The best current choice for support-wall material is graphite, which is a solid up to 
about 7200° R (4000° K). Thus, for a practical system, the wall temperature must be 
maintained below 7200° R (4000° K) with acceptable wall cooling ratios. Figure 6 indi- 
cates that at the reference value of 77 = 0. 25 the support wall must operate at 9830° R 
(5460° K), and that this temperature cannot be brought into the range of realistic material 
limits even by increasing the wall cooling fraction to 77 = 0. 4 at the loss of rocket spe- 
cific impulse (see table I). It is thus clear that the feasibility of this radiating liquid-fuel 
rocket depends on a scheme for reducing the support-wall temperature at acceptable wall 
cooling rates. 


Fuel-Mixture Temperatures With Zoned Heat Sources And 
Time Dependence Of Zoned Fuel Concentrations 

Within the restrictions imposed by system requirements, the wall temperature may 
be lowered considerably by zoning the fissioning fuel into a region in the liquid near the 
gas interface. Without considering for the moment the details of confining the fuel to a 
small region in the liquid, the steady-state temperatures with a step distribution of heat 
sources or fuel concentration within the liquid can be calculated from equations (19), (20), 
and (21). Figure 7 was constructed by using the calculated turbulence coefficient A = 

11. 0 and a typical zone fraction f3 = 0. 90 (for fuel confined to one-tenth of the liquid 
thickness). A comparison of figures 6 and 7 indicates that, by fuel zoning with (3 = 0.90 
(and no fuel diffusion), support-wall temperatures are reduced about 1200° R (670° K) 
below those without zoning. Making the fuel region narrower, such that |3 > 0. 90, would 
not measurably enhance these results. 

Since fuel zoning is a potential means of reducing the wall temperature, the ability 
to actually confine nuclear fuel to a small fraction of the liquid thickness must be inves- 
tigated. The fuel is initially solid and is zoned. Upon liquidification of the fuel mixture 
at rocket startup, steep radial fuel concentration gradients cause diffusion of fuel away 
from the zoning region and result in redistribution of fuel throughout the liquid mixture. 
The rate of this fuel redistribution is evaluated using the analysis of the TIME DEPEN- 
DENCE OF ZONED FUEL CONCENTRATIONS section. 


15 



■*■■ II ■■■ Mil 


The time dependence of fuel concentration with initial step distribution is calculated 
using equations (27) and (29). Normalized concentration profiles for the typical zoning 
fraction, /3 = 0. 90, are shown for several dimensionless times after rocket startup in 
figure 8. Of particular interest is the time fljyg wlien the fuel concentration at the wall 
reaches one -half the equilibrium value, that is, when C Q / C = 0. 5. 

The quantity 0jy ^ is assumed to be a measure of the time after which the advan- 
tages of reduced wall temperature due to fuel zoning are lost. Figure 9 shows the time 
variation of fuel concentration at the solid wall for various /9 values. Values of fuel 
concentration were calculated using equation (29). The maximum duration of zoning 
effects is taken from figure 9 to be the largest possible value of e( l ua l to 0. 14, 

which corresponds to the case where all the fuel heat sources are initially zoned into a 
very thin layer at the gas interface, or j8 s 1. 0. With known mass diffusion coefficient 
this limit corresponds to tjyg = (L^fljyg/D) = 0. 14 (L^/D). 

Assuming only molecular mass diffusion occurs and using the estimated value of D 
given in the appendix give the maximum duration of zoning effects as tjyg = 25 minutes. 
In reality, turbulent effects, which were shown previously to be an important mechanism 
of heat transfer, should be included in estimating the value of D. Thus when turbulence 
is accounted for, the maximum expected value of t^yg becomes so small that any ad- 
vantages of wall temperature reduction due to initial fuel zoning are expected to be lost 
within a fraction of realistic rocket operating life. On the basis of a simple mass dif- 
fusion model, it is thus concluded that initial fuel zoning cannot be expected to lower 
significantly the solid wall temperature over rocket operating life. 

Several means of inhibiting the rapid redistribution of zoned fuel may be considered. 
If the NbC solvent, which has a lower density than the UCg fuel, is replaced with a 
material having the same physical properties as NbC but a greater density than the UCg 
fuel, it is conceivably possible that the rotating fuel element would act like a centrifuge. 

In effect, the lower density fuel would tend to collect in the region near the gas -liquid 
interface as desired. 

An estimate of this effect was made using the centrifuge equation given in Bird, et al. 
(ref. 6). For the design radial acceleration of 387 g Q , a fuel mixture depth of 0. 0095 
meter, and the extreme case of a solvent having twice the density of UCg, the equilibrium 
fuel concentration in the centrifuge case was found to be only 6 percent greater at the gas 
interface than at the solid wall. Thus the centrifuge effect would not significantly inhibit 
the rates at which zoned fuel is redistributed by mass diffusion. 

Another possible means of confining the fuel would be to float it on a denser sub- 
stance which, although liquid, would not be wetted by the liquid fuel. Thus, surface ten- 
sion effects would tend to confine the fuel to a layer resting on another liquid layer. How- 
ever, even if such a more dense nonwetting substance could be found for this application, 
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turbulence, which would be induced in the liquid fuel and its support substance, would 
probably cause breakup and mixing of the fueled and unfueled layers, again redistributing 
the fuel and heat sources throughout the liquid. 

Thus, for the possibilities just considered, it is clear that with initial fuel zoning the 
fuel redistributes itself throughout the liquid region of the fuel element shortly after 
rocket startup. Therefore, for the greater part of rocket operating life the fuel temper- 
ature distribution would be as shown in figure 5 . 

Support-Wall Temperature With an Insulating Vapor 
Layer at the Liquid-Solid Interface 

An alternative means of reducing support-wall temperature is by means of a vapor 

layer built up at the liquid-solid interface. The analytic details of generating such a 

vapor layer, whether by sublimation of a graphite support wall, or injection through a 

porous wall, are not treated here. However, a simple derivation in the analysis section 

leading to equation (31) indicates that the maximum thickness of a stable vapor layer that 

\ 

remains attached to the support wall is about 0. 57 millimeter. 

The temperature drop across a vapor layer was evaluated using the radiative heat- 
transfer analysis given in reference 4. The support -wall temperature was calculated for 
vapor layers of various optical thicknesses by using the values of liquid -mixture temper- 
ature T q and heat flux q” at the wall taken from figure 5. The results are shown in 
figure 10. The actual value of optical thickness x corresponding to the vapor -layer 
thickness of 0. 57 millimeter is difficult to estimate as pointed out in the analysis; how- 
ever, a representative value was calculated assuming a carbon vapor layer at 9000° R 
(5000° K). The calculated value of x = 1. 8 is possibly an overestimate and should be 
experimentally verified. Support-wall temperatures for the assumed vapor layer are 
given by the curve for x = 1. 8 in figure 10. These estimates for a stable vapor layer 
at the support wall indicate a wall temperature of 8500° R (4720° K) for the design value 
77 = 0. 25. Also, a wall temperature of 7200° R (4000° K) may be achieved by increasing 
wall cooling to 77 = 0. 32, which decreases specific impulse to 1250 seconds. 


Engine Performance Feasibility 

The high support-wall cooling rates require the propellant to be preheated to some 
inlet temperature before entering the fuel element proper. Reference 1 gives the rela- 
tion that results between rocket specific impulse and wall cooling ratio for the case where 
propellant inlet temperature is limited to a 3500° R (1940° K) maximum by rocket 
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structural material limits. This relation is given in table I. 

The support-wall temperatures for the three principal cases considered are replotted 
from previous figures as functions of wall cooling ratio and the corresponding specific 
impulse in figure 11. The curve for fuel zoning is mainly of theoretical interest and can- 
not be physically realized due to rapid fuel redistribution rates. 

Two points on these curves are especially important. For the case of unzoned fuel 
directly on the support wall, figure 11 indicates that, for the reference specific impulse 
of 1430 seconds, solid support-wall materials with structural limits above 9830° R 
(5460° K) are required. The curve for the insulating vapor layer, while requiring fur- 
ther experimental verification, indicates that current wall materials, with a 7200° R 
(4000° K) limit, may permit operation with a reduced specific impulse of 1250 seconds. 


CONCLUSIONS 

This analysis of heat transfer in the liquid fuel mixture extends the work of refer- 
ence 1 on a new radiating liquid-core nuclear rocket concept, which consists of a liquid 
carbide fuel mixture supported on the inside wall of a rotating tube. The axially flowing 
hydrogen propellant is heated by thermal radiation from the liquid surface. An expres- 
sion for combined turbulent and molecular heat transfer in the liquid fuel was developed 
and then used to calculate temperatures at the solid support wall for a 1430-second spe- 
cific impulse rocket with 9500° R (5280° K) fuel surface temperature. 

With turbulent mixing in the liquid fuel accounted for, support -wall temperature for 
the 1430-second rocket was 9830° R (5460° K), well in excess of practical wall material 
limits. 

With fuel zoned into one-tenth of the liquid thickness farthest from the support wall, 
it was found that wall temperatures are reduced about 1200° R (660° K) by fuel zoning. 
However, a mass diffusion analysis indicates that fuel zoning could not be maintained 
over rocket lifetime. Thus, the advantages of reduced wall temperature due to fuel zon- 
ing cannot be physically realized. 

The insulating effect of a stable vapor layer built up at the solid-liquid interface was 
also evaluated. The estimated value of radiation optical thickness for such a vapor layer 
is probably not very accurate; however, for the particular case studied it was found that 
wall temperatures were significantly reduced by a vapor layer at the wall. A specific 
impulse of 1430 seconds required a 8500° R (4720° K) support -wall temperature: a spe- 
cific impulse reduced to about 1250 seconds would require a more realistic 7200° R 
(4000° K) wall temperature. 

The final conclusions reached are as follows: 

1. Turbulent mixing in the liquid fuel must be accounted for in the estimates of 
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support -wall temperature. 

2. To realize a design specific impulse of about 1430 seconds, solid support-wall 
materials with melting points near 9830° R (5460° K) must be developed. 

3. The buildup of a vapor layer at the wall may permit use of current materials, 
such as graphite, for a support wall at 7200° R (4000° K), but would probably reduce the 
attainable specific impulse to about 1250 seconds. The acutal hydrodynamic and radiative 
heat-transfer characteristics of such an insulating vapor layer require further experi- 
mental verification. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 12, 1967, 

122-28-02-18-22. 
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APPENDIX 

LIQUID-FUEL-MIXTURE PROPERTIES 

For determination of properties, the liquid fuel mixture is assumed to consist of 
liquid NbC at 9000° R (5000° K). Property data were taken from the following sources, 
and in most cases were estimated using semi-empirical relations from the indicated 
sources: 

. 6750 (3750) 

. 8570 (4770) 

13 000 (7200) 
.... 7150 
585 

- 4. 75xl0~ 3 

.... 43.3 
. 6.77X10" 7 
. . . 0.0647 

.... 1.49 
. 8. 37xlO' 9 


Melting point temperature, °R (°K) (ref. 7) 

Boiling point temperature (for 1 atm (1. 01x10 N/m )), 
°R (°K) (ref . 7) 


Boiling point temperature (for 200 atm (2. 02xl0 7 N/m^)), 


°R (°K) (ref. 8) 

Density (estimated using ref . 7), p^, kg/m 

Specific heat (estimated using ref. 9), , J/(kg)(°K) . • . . 

L 2 

Viscosity (estimated using ref. 10), ju^, (N)(sec)/m 

Molecular thermal conductivity (estimated using ref. 9), k^, 

J/(m)(sec)(°K) 

Kinematic viscosity (t> L = (p L /p L )), m /sec 

Prandtl number (Pr^ = (p^Cp ^ y. Pr^ 

L l 7 j 

Surface tension (estimated using ref . 11), a, N/m 

Molecular mass diffusion coefficient (estimated using ref. 12), 
D, m /sec 
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(a) Conceptual representation of radiating fuel element. 
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Dimensionless distance from 
support wall, y/L 

Figure 2. - Variation of molecular plus turbu- 
lent thermal conductivity with distance from 
support wall. Ratio of turbulent to molecular 
thermal conductivity, A = 11,0. 


(b) Nomenclature and general characteristics of analyzed fuel element. 
Figure 1. - Liquid-core fuel element of nuclear rocket. 
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Figure 3. - Schematic of vapor bubble on support 
wall showing nomenclature. 
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Figure 4. - Temperature profile in liquid fuel mixture for 
various wall cooling ratios with molecular conduction 
only and uniform heat source distribution. Ratio of 
turbulent to molecular thermal conductivity, A = 0; 
fraction of liquid thickness which has no heat sources 
(3 = 0. 
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Figure 5. - Temperature profile in liquid fuel mixture for 
various wall cooling ratios with molecular plus turbu- 
lent heat transport and uniform heat source distribu- 
tion. Ratio of turbulent to molecular thermal conduc- 
tivity, A = 11. 0; fraction of liquid thickness which has no 
heat sources, (3 = 0. 




Wall cooling ratio, r\ 


Figure 6. - Variation of support -wall temperature with 
wall cooling ratio for various turbulence coefficients 
and uniform heat source distribution. Fraction of 
liquid thickness which had no heat sources, (3 = 0. 
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Figure 7. -Temperature profile in liquid fuel mixture for 
various wall cooling ratios with molecular plus turbulent 
heat transport and zoned heat sources. Ratio of turbulent 
to molecular thermal conductivity, A = 11.0, fraction of 
liquid thickness which had no heat sources, 3 = 0. 9. 


Dimension- 
less time, 

0 = Dt/L 2 



Dimensionless distance from support wall, y/L 

Figure 8. - Fuel concentration variation with dimensionless time 
with fuel zoned initially into one-tenth of liquid thickness. 
Fraction of liquid thickness which had no heat source, {3 = 0. 9. 
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Figure 9. - Variation of fuel concentration at support wall 
with time for various zone fractions. 
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Figure 10. - Variation of support-wall temperature with 
wall cooling ratio for vapor layers of various optical 
thicknesses at wall. Ratio of turbulent to molecular 
thermal conductivity, A = 11.0; fraction of liquid 
thickness which had no heat sources, (3 = 0. 
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Specific impulse, sec 

Figure 11. - Variation of support-wall temperature with 
wall cooling ratio and specific impulse for various cases 
analyzed and turbulent mixing included. 
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